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NOMENCLATURE 


A 

Surface  area  of  sphere    (ft  ) 

a 

Amplitude  of  oscillation  (ft) 

P 

Specific  heat  (Btu/lb>m-F) 

D 

Diameter  of  sphere  (ft) 

f 

Frequency  of  oscillation  (cyc/hr) 

h 

Mean  unit  surface  conductance  (Btu/hr-ft 

2-F) 

"hr 
*  c 

Mean  unit  surface  convective  conductance 

(Btu/hr-f t2-F) 

"h- 
r 

Mean  unit  surface  radiative  conductance 

(Btu/hr-ft^-F) 

k 

Thermal  conductivity    (Btu/hr-f t-F) 

m 

Mass  of  the  sphere  (lt>ra) 

T 

Absolute  temperature  (R) 

Tm 

Mean  temperature  (R) 

L 

Film  temperature    (Tm  +  ToJ/2  (R) 

To 

Temperature  when  8  ■  0  (R) 

Tfl 

Temperature  when  8  ■  0C  (R) 

T 

Temperature  of  ambient  fluid  (R) 

V 

Volume  of  sphere    (ft  ) 

6 

Average  emissivity  of  surface 

Time  (hr) 

Dvnamic  viscositv     ^lh  /ft-hr} 
■* j itoiHi^  v  locuo i  y     *    m        u & j 

r 

Density  (lb™/ft^ 

S  f"  ^  nH  2*  n  -Rn  1  f  yma  n    <-»  nn  c  t-  o  n  t-       /A    171  A    v    1 A*^ 

•Jtc^iidii  Duit2iuan  constant     ^u.i/m  x  1U 

otu/nr-rtfc-R  ) 

a) 

Angular  velocity  (1/hr) 

Bi 

Biot  number  (hV/kA) 

Nu 

Nusselt  number  (hcD/kf) 

Re 

Reynolds  number    (  p  Daio/^) 

INTRODUCTION 


In  recent  years,  much  work  has  been  done  with  regard  to  the  heat  trans- 
fer from  spheres.    The  investigations  thus  far  have  been  concerned  with 
natural  convection,  combined  natural  and  forced  convection  from  stationary 
spheres,  and  convection  from  a  rotating  sphere.    This  report  presents  the 
results  of  an  experimental  study  of  the  heat  transfer  from  an  oscillating 
sphere  to  a  quasi  infinite  medium  as  a  function  of  frequency  and  amplitude. 

EXPERIMENTAL  EQUIPMENT 

Plate  I  is  a  photograph  of  the  equipment  used  in  the  experiments.  The 
one  inch  diameter  sphere  was  made  of  copper  and  was  within  0.005  inches  of 
being  spherical.    The  size  and  material  were  selected  so  that  the  internal 
thermal  resistance  was  negligible  compared  to  the  thermal  resistance  at 
the  surface  of  the  sphere.    Kreith  (A)  states  that  when  the  internal  resis- 
tance is  less  than  ten  percent  of  the  external  resistance,  the  error  in 
assuming  the  temperature  throughout  the  sphere  to  be  constant  is  less  than 
five  percent.    That  is,  Bi  <Co.l. 

The  mechanism  used  to  obtain  sinusoidal  oscillations  of  the  sphere  in 
a  vertical  direction  was  a  Scotch  yoke.    The  drive  system  consisted  of  a 
1/15  horsepower  variable  speed  a-c  motor  with  a  speed  reducing  pulley.  A 
variable  length  crank  arm  connected  the  pulley  shaft  to  the  Scotch  yoke. 
The  frequency  of  oscillation  was  varied  by  using  an  auto transformer.  The 
low  torque  characteristic  of  the  motor  at  low  input  voltage  made  it  necessary 
to  balance  the  system  in  order  to  obtain  true  sinusoidal  oscillations.  At 
frequencies  less  than  three  cps,  the  frequency  varied  during  a  test  as  much 
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as  four  percent  from  a  mean  value.    Above  three  cps,  a  stroboscope  showed 
that  the  frequency  remained  constant  over  the  entire  time  of  the  test  run. 

The  sphere  was  connected  to  the  Scotch  yoke  by  means  of  a  six  inch 
length  of  1/8  inch  diameter  Megopak  thermocouple  wire  and  a  3/8  inch  dia- 
meter thin  walled  steel  tube.    The  Megopak  was  soldered  to  the  sphere  with 
1200  F  silver  solder  so  that  the  iron-cons tantine  thermocouple  junction 
was  located  at  the  center  of  the  sphere.    The  sphere  and  drive  system 
could  be  lowered  from  the  position  shown  in  plate  I  for  immersion  in  a 
water  filled  barrel  to  a  depth  of  approximately  six  inches. 

The  sphere  was  heated  by  two  resistance  heaters  located  on  opposite 
sides  of  a  nine  inch  cube  shaped  oven.     The  oven  was  made  in  halves  and 
spring  loaded  so  that  it  would  open  quickly  from  around  the  sphere  and 
swing  away  a  distance  of  approximately  three  feet  (see  plate  II). 

A  Honeywell  recorder  was  used  to  obtain  the  temperature -time  records 
for  the  sphere  cooling  in  air.  When  the  cooling  medium  was  water,  a  San- 
born d-c  amplifier  and  recorder  was  used  in  order  to  obtain  a  chart  speed 
of  100  mm  per  second.  A  disadvantage  of  this  instrument  was  that  its  ac- 
curacy is  not  guaranteed  for  input  voltages  less  than  5  mv. 

EXPERIMENTAL  TECHNIQUE 

Each  test  was  begun  by  heating  the  sphere  to  a  temperature  well  above 
that  desired  for  the  upper  limit  of  the  cooling  range.    During  heating, 
the  sphere  was  brought  to  the  desired  frequency  of  oscillation.    Then  the 
oven  was  triggered  and  a  temperature -time  history  for  the  sphere  was  re- 
corded. 
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For  the  tests  in  air,  the  oven  was  closed  immediately  in  order  to  pre- 
vent radiation  from  the  oven  to  the  sphere  and  heating  of  the  air.    This  was 
done  while  the  sphere  temperature  was  still  well  above  the  temperature  cor- 
responding to  time  zero  for  the  test.    Thus,  the  cooling  process  was  not 
affected  by  the  oven  during  the  test  period.    For  the  tests  in  water,  the 
heating  to  a  higher  temperature  eliminated  the  effects  of  dunking. 

The  cooling  of  a  sphere  with  negligible  internal  resistance  is  analogous 
to  the  discharge  of  a  capacitor  in  an  R-C  electrical  circuit  (A) .    A  first 
law  analysis  on  the  sphere  yields  the  equation 

(1)  -  mepdT  -  hA(T  -  T60)d9 

The  solution  for  this  differential  equation  gives  the  temperature  of  the 
sphere  as  a  function  of  time.    It  is 

(2)  TQ  -  Tw  ^ 

mc« 

The  term  «—*  is  called  the  time  constant  and  will  be  designated  by  6C  in 

hA 

this  report.  When 

mcp 

(3)  9  «=  9C  -  — 

hA 

equation  (2)  reduces  to 

<*>         Tqc  -  T„-  l/e(  TQ  -  T^  ) 
The  temperature  range  for  each  test  was  determined  from  equation  (4) 
by  choosing  an  initial  temperature  TQ  and  calculating  the  final  temperature 
Tec-    Because  the  room  temperature  T    changed  from  day  to  day,  initial 
temperatures  of  610  F,  590  F,  and  580  F  were  selected  for  the  tests  in  air 
in  order  to  obtain  a  temperature  range  less  than  that  of  the  recorder.  An 
initial  temperature  of  796  F  was  selected  for  the  tests  in  water  to  include 
the  four  regimes  of  cooling  on  the  temperature-time  record.    The  time 
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constant  for  each  test  run  was  calculated  from  the  temperature -time  record 
by  dividing  the  measured  distance  between  TQ  and  Tqc  by  the  chart  speed  of 
the  recorder.    With  the  time  constant  known,  the  unit  surface  conductance 
was  calculated  from  equation  (3) . 

Since  the  specific  heat  of  copper  varies  with  temperature  (B),  a  value 
corresponding  to  the  mean  temperature  of  the  sphere  over  the  test  range  was 
used  in  equation  (3).    The  mean  temperature  was  determined  from  the  equation 

re  -jjj 

/(To  -  T„  )Q  de 
9 

When  equation  (5)  is  integrated  from  0  to  9c,  it  reduces  to 

<6>         Tm  "  T«  =  (1  "  1/«XT0  " 
It  was  noted  that  the  value  of  h  obtained  by  using  a  value  of  cp  corresponding 
to  an  average  temperature  calculated  by  the  equation 

(7)         T        -    T°  +  T*» 
V'         Aavg.  "   5  


was  less  than  three  percent  lower  than  the  "h  obtained  using  cp  corresponding 
to  Tm. 

The  heat  loss  to  the  Megopak  by  conduction  was  determined  to  be  less  than 
one  percent  by  calculation  of  both  sides  of  equation  (1)  in  the  form 

(8)  mcp(T0  -  T9c)  -  hAO^  -  ?„)9C 

The  heat  transfer  by  radiation  was  found  to  be  a  major  part  of  the 
total  heat  transfer  for  the  tests  in  air;  however,  it  was  negligible  when 
the  cooling  medium  was  water.    The  unit  surface  conductance  was  separated 
into  two  parts  by  the  equation 

(9)  h  =  hr  +  hc 
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The  unit  radiative  conductance  was  calculated  from  the  equation 

_     €<r  at  -  it  > 

(10)  h  =   

(Tm  -  T*  > 

and  substituted  into  equation  (9)  to  determine  hc  for  each  test  run.  Equation 
(10)  is  of  the  form  given  by  Yuge  (D) .    The  emissivity  of  the  copper  was 
taken  as  0.87  from  Kreith  (A). 

RESULTS 

Several  tests  were  made  in  air  at  zero  frequency.    The  values  of  hc 
obtained  agreed  with  the  correlation  equation  published  by  Nordlie  and 
Kreith  (C)  within  their  prescribed  limits.    This  check  served  to  establish 
the  reliability  of  the  test  equipment  and  calculation  procedure. 

In  plate  HI,  the  cooling  curve  for  the  sphere  oscillating  at  8  cps 
and  one  inch  amplitude  is  compared  to  that  for  the  sphere  at  zero  fre- 
quency.   Also  shown  is  a  plot  of  equation  (2)  using  the  h  obtained  from 
test.    These  curves  are  typical  of  those  for  the  series  of  tests  on  the 
sphere  cooling  in  air.    The  series  covered  a  range  of  frequencies  from  0 
to  10  cps. at  amplitudes  of  1,  3/4,  and  1/2  inch.  A  plot  of  h  versus  f  is 
shown  in  plate  IV.    The  mathematical  relation  for  h  as  a  function  of  f 
in  cps  for  each  amplitude  was  determined  by  the  calculus  of  finite  dif- 
ferences.   The  relations 

h  -  3.66  +  0.175  f  +  0.015  f2         (a  -  1  inch) 

(11)  h  -  3.66  +  0.0655  f  +  0.0145  f2      (a  -  3/4  inch) 
h  -  3.66  +  0.045  f  +  0.005  f2         (a  -  1/2  inch) 

describe  the  data  to  within  ±  5  percent.    A  plot  of  Nusselt  number  versus 
Reynolds  number  for  the  data  obtained  is  shown  in  plate  V.    The  fluid  pro- 
perties were  evaluated  at  the  film  temperature. 
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Although  no  equation  is  given,  the  curve  shown  correlates  the  data  to  with- 
in +  15  percent . 

Four  tests  were  made  at  an  amplitude  of  1/2  inch  using  water  at  60  F 
as  the  cooling  medium.    The  cooling  curves  for  these  tests  are  shown  in 
plate  VI.    In  each  test,  the  sphere  was  cooled  from  796  F  to  approximately 
80  F.    Four  regimes  of  cooling  are  covered  by  these  curves  (B) .    The  re- 
gimes are  indicated  approximately  for  the  curve  with  a  frequency  of  9.2  cps. 
During  film  boiling,  the  plot  shows  that  equation  (2)  is  correct  to  within 
5  percent  and  that  the  unit  surface  conductance  is  not  a  function  of  fre- 
quency.   However,  it  was  noted  that  the  transition  from  film  boiling  to 
nucleate  boiling  occured  at  a  higher  temperature  for  a  higher  frequency 
of  oscillation.    A  visual  observation  indicated  that  the  sphere  moved  out 
of  the  film. 
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The  heat  transfer  from  a  one  inch  diameter  copper  sphere  oscillating  in  a 
quasi  infinite  medium,  both  air  and  water,  has  been  studied.    Unit  surface  con- 
ductances were  experimentally  determined  for  frequencies  from  0  to  10  cps  at 
amplitudes  of  1,  3/4  and  1/2  inch.    The  oscillating  conditions  correspond  to 
a  Reynolds  number  range  from  33  to  1.6  x  10^.    The  results  of  the  tests  in  air 
are  shown  graphically  by  a  plot  of  Nusselt  number  versus  Reynolds  number.  The 
curve  shown  correlates  the  results  to  within  +  15  percent.    Four  tests  were 
run  in  water  using  an  amplitude  of  1/2  inch  and  frequencies  of  0,  2.5,  6.33, 
and  9.2  cps.    The  cooling  curves  show  that  the  unit  surface  conductance  is 
not  a  function  of  frequency  during  film  boiling  and  that  the  transition 
from  film  boiling  to  nucleate  boiling  occurs  at  a  higher  temperature  for  a 
higher  frequency  of  oscillation. 


